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Solar neutrino measurements from 1258 days of data from the Super-Kamiokande detector are 
presented[? ]. The measurements are based on recoil electrons in the energy range 5.0-20.0 MeV. 
The measured solar neutrino flux is 2.32 ± 0.03 (stat.) 1° q? ( s y s -) x 10 6 cm -2 s _1 , which is 
45.1 ± 0.5 (stat.) ±{\l (sys.)% of that predicted by the BP2000 SSM. The day vs night flux asym- 
metry ($„ — <&d)/<&average is 0.033 ± 0.022 (stat.) 1q qi2 ( s y s 0- The recoil electron energy spectrum 
is consistent with no spectral distortion (x 2 /d.o.f. = 19.0/18). The seasonal variation of the flux 
is consistent with that expected from the eccentricity of the Earth's orbit (x 2 /d.o.f. — 3.7/7). For 
the hep neutrino flux, we set a 90% C.L. upper limit of 40 x 10 3 cm -2 s _1 , which is 4.3 times the 
BP2000 SSM prediction. 



This preprint is almost identical to the report submitted to Phys- 



ical Review Letter. We have added to this preprint a few tables of 
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Solar neutrinos have been detected using chlorine-, 
gallium-, and water-based detectors jl], g, ||, ||, |J; all 
have measured significantly lower solar neutrino flux than 
predicted by Standard Solar Models (SSMs) @, @, §• 
This disagreement between the measured and expected 
solar neutrino flux, known as the "solar neutrino prob- 
lem" , is generally believed to be due to neutrino flavor 
oscillations. Signatures of neutrino oscillations in Super- 
Kamiokande (SK) might include distortion of the recoil 
electron energy (E rccoi i) spectrum, difference between the 
night-time solar neutrino flux relative to the day-time 
flux, or a seasonal variation in the neutrino flux. Obser- 
vation of these effects would be strong evidence in sup- 
port of solar neutrino oscillations independent of abso- 
lute flux calculations. Conversely, non-observation would 
constrain oscillation solutions to the solar neutrino prob- 
lem. We describe here solar neutrino measurements from 
1258 days of SK data. 

SK, located at Kamioka Observatory, Institute for Cos- 
mic Ray Research, University of Tokyo, is a 22.5 kton 
fiducial volume water Cherenkov detector that detects 
solar neutrinos via the elastic scattering of neutrinos off 
atomic electrons. The scattered recoil electron is de- 
tected via Cherenkov light production, allowing both the 
direction and total energy to be measured. These quan- 
tities are related to the original neutrino direction and 
energy. Detailed descriptions of SK can be found else- 
where § | 

The 1258-day solar neutrino data were collected in four 
periods with different trigger thresholds between May 
31, 1996 and October 6, 2000 (table |). The analysis 
threshold has been at 5.0 MeV except for the first 280 
days where the data were analyzed with a threshold of 
6.5 MeV. The analysis threshold is determined by the 
level of irreducible background events and the event trig- 
ger threshold. An event is triggered when the sum of 
PMTs registering a hit in a 200 nsec time window (Nhu) 
is above a threshold (table |) . This threshold should be 
sufficiently low that the trigger efficiency at the analy- 
sis threshold is nearly 100%. The lowering of the trigger 
threshold in periods 2-4 was made possible by the ad- 
dition of a software filter to the data acquisition system 
that removes a large portion of background events. This 
removal is accomplished by reconstructing the event ver- 
tex and rejecting events with vertices within 2 m of the 
inner detector wall, most of which are due to external ra- 
dioactivity. Each lowering of the trigger threshold in the 
course of the experiment was made possible by increasing 
the number of computers that run the filter program. 

There are 2.0 x 10 9 events in the raw data sample before 
background reduction. After removing cosmic ray muon 
events, the sample in the 22.5 kton fiducial volume with 



numerical values that were omitted from the PRL draft. 



Run 


N hit 


50/95% 


Analysis 


Live-time 


period 


threshold efficiency threshold 


(days) 






(MeV) 


(MeV) 




(1) May 1996 <- 


- 40.6 


5.7 / 6.2 


6.5 


280 


(2) May 1997 <- 


- 34.5 


4.7 / 5.2 


5.0 


650 


(3) Sep. 1999 - 


- 30.4 


4.2 / 4.6 


5.0 


320 


(4) Sep. 2000 - 


- 27.7 


3.7 / 4.2 


5.0 


8 



TABLE I: The trigger and analysis thresholds and live-times 
during which they were used. The third column shows the 
recoil electron energy at which the trigger is 50% and 95% 
efficient. The software filter was added starting in May 1997. 



energy between 5.0-20.0 MeV contains 3.0 x 10 7 events. 
The dominant background sources in the low-energy re- 
gion (E < 6.5 MeV) are 222 Rn in the water and external 
radioactivity; in the high-energy region (E > 6.5 MeV), 
radioactive decay of muon-induced spallation products 
accounts for most of the background. Background reduc- 
tion takes place in the following steps: first reduction, 
spallation cut, second reduction, and external gamma- 
ray cut. The first reduction includes cuts that remove 
events due to electronic noise and arcing PMTs. In ad- 
dition, a cut on the goodness of the reconstructed vertex 
is used to remove obvious background events originating 
from various non-physical sources. The number of re- 
maining events after the first reduction is 1.5 x 10 7 . The 
spallation cut has been improved compared to that used 
in earlier publications [|[ [To], |ll| . We have improved the 
likelihood functions used in removing spallation events 
and introduce a new cut for 16 N events that originate 
from absorption of cosmic ray stopped [iT on 16 O. The 
number of events in the high-energy region (6.5-20 MeV) 
before and after spallation cut is 1.6 x 10 6 and 3.3 x 10 5 , 
respectively. The spallation cut is 79% efficient for solar 
neutrino events. The second reduction removes events 
with poor vertex fit quality or with blurred Cherenkov 
ring patterns, characteristics of low-energy background 
events and external gamma rays. This newly introduced 
reduction step has improved the signal-to-noise ratio in 
the low-energy region by almost an order of magnitude. 
The number of events before and after the second reduc- 
tion in the 5.0-6.5 MeV region are 1.0 x 10 7 and 1.4 x 10 6 
events, respectively. In addition, the gamma-ray cut, 
which removes external events, has been tightened for 
those events with E < 6.5 MeV. The combined efficiency 
of the first reduction, second reduction, and the exter- 
nal gamma ray cut for solar neutrino events is ~ 73% 
for E > 6.5 MeV, and - 52% for E < 6.5 MeV. After 
these reduction steps, 236,140 events remain in the fidu- 
cial volume above 5 MeV, with S/N « 1 in the solar 
direction. 

The SK detector simulation is based on GEANT 
3.21 [1^]. The energy scale was measured using a larger 
sample of data from an in situ electron linear accelera- 
tor |)| (LINAC) compared to that used in earlier results. 
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The detector simulation's reliability was tested using the 
well-known f3 decay of 16 N, which is produced in situ by 
an (n,p) reaction on 16 0. Fast neutrons for this reaction 
are produced using a portable deuterium-tritium neutron 
generator (DTG) jl3| . The energy scale measured by the 
DTG agrees with that from the LINAC within ±0.3%. 
The total systematic uncertainty in the absolute energy 
scale, including possible long term variation and direction 
dependence, is ±0.6%. 

We compare our solar neutrino measurements against 
reference fluxes and neutrino spectra in order to search 
for signatures of neutrino oscillations. For E Teco u > 
5.0 MeV, solar neutrinos are expected to come almost 
exclusively from the decay of 8 B, with a slight ad- 
mixture of neutrinos from 3 He-proton (hep) fusion. For 
the absolute flux of 8 B and hep neutrinos, we take the 
BP2000 H SSM as our reference. The /3 decay spectrum 
of the 8 B neutrinos is dominated by the transition to a 
broad excited state of 8 Be, which decays immediately to 
two a particles. Bahcall et al. |l4j use a neutrino spec- 
trum deduced from a comparison of world data on 8 Be 
a decay [[IBJ |l^, 17 with the direct measurement of the 
positron spectrum from 8 B decay measured by Napoli- 
tano, Freedman, and Camp fl8|| . Energy-dependent sys- 
tematic errors are deduced from a combination of exper- 
imental uncertainties and the theoretical uncertainties in 
radiative and other corrections that must be made to 
convert the charged particle data into a neutrino spec- 
trum JlJ]. Recently, Ortiz et al. Jl§ have made an im- 
proved measurement of the 8 B spectrum based on 8 Be a 
decay in which some of the major sources of systematic 
errors present in previous measurements were reduced 
or eliminated. We have adopted the neutrino spectral 
shape and experimental uncertainties from this measure- 
ment. These experimental uncertainties were then added 
in quadrature with the theoretical uncertainties given by 
Bahcall et al. Q. 

The solar neutrino signal is extracted from the data 
using the cosf? sun distribution (Fig. |l|). The angle 9 sun is 
that between the recoil electron momentum and the vec- 
tor from the sun to the Earth. The solar neutrino flux is 
obtained by a likelihood fit of the signal and background 
shapes to the cos 9 sun distribution in data. The signal 
shape is obtained from the known angular distribution 
and detector simulation, while the background shape is 
nearly flat in cos6> sun . In the 8 B flux measurement, the 
data are subdivided into 19 energy bins in the range 5.0— 
20.0 MeV (binning as in Fig. |). The likelihood function 
is defined as follows: 



19 - ( Y r S +Bj ) 



n 



N 3 \ 



Yl [Bj ■ F b (cos 6i ,Ei) + Y 3 -S- F s (cos ft , E t 



(1) 

S is the total number of signal events, while Nj, Bj, and 
Yj represent the number of observed events, the number 
of background events, and the expected fraction of signal 
events in the j-th bin, respectively. Ff, and F s are the 
probability for the background and signal events as a 
function of cos0 sun and energy (Ej) of each event. The 
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FIG. 1: cos6> SUI1 distribution for reconstructed energy 
E = 5.0-20.0 MeV. The points represent observed data, the 
histogram shows the best-fit signal level plus background, and 
the nearly horizontal line shows the estimated background. 
The peak at cos 8 sun = 1 is due to solar neutrinos. 



likelihood function is maximized with respect to S and 
Bj. For the energy spectrum measurement, each term in 
the product over bins is maximized separately. 

The best-fit value of S is 18, 464±204 (stat.)t^ (sys.), 
which is 45.1 ± 0.5 (stat.) (sys.)% of the reference 
prediction. The corresponding 8 B flux is: 



2.32 ± 0.03 (stat.)lo:o7 (sys.) 



x 10 b 



The total systematic error is _3o%> with the largest 
sources coming from the reduction cut efficiency (^'Zy), 
energy scale and resolution (±1.4%), systematic shifts in 
the event vertex (±1.3%), and the angular resolution of 
the recoil electron momentum (±1.2%). 

Fig. |^ shows the solar neutrino flux as a function of 
the solar zenith angle 9 Z (the angle between the vertical 
axis at SK and the vector from the sun to the Earth). 
Numerical values are shown in Table |l|. The day-time 
solar neutrino flux <&d is defined as the flux of events 
when cos 9 Z < 0, while the night-time flux <!>„ is that 
when cos 9 Z > 0. The measured fluxes are: 

$ d = 2.28 ±0.04 (stat.)lg;o? MO x 1()6 cm-V 1 
= 2.36 ±0.04 (stat.)l°;° 8 (sys.) x 10 6 cm^s" 1 

Some neutrino oscillation parameters predict a non-zero 
difference between and $^ due to the matter effect in 
The Earth's mantle and core (2(J. The degree of this dif- 
ference is measured by the day-night asymmetry, defined 

as .4= ($n-$d)/®average, where ^average = |(*n + *d)- 

We find: 



A = 0.033 ± 0.022 (stat.)!";^ 3 (sys.) 

Including systematic errors, this is 1.3 a from zero asym- 
metry. Many sources of systematic errors cancel out 
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FIG. 2: The solar zenith angle (t9 z ) dependence of the solar 
neutrino flux (error bars show statistical error). The width 
of the night-time bins was chosen to separate solar neutrinos 
that pass through the Earth's dense core (cos 9 Z > 0.84) from 
those that pass through the mantle (0 < cos(9 z < 0.84). The 
horizontal line shows the flux for all data. 
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FIG. 3: The measured 8 B + hep solar neutrino spectrum 
relative to that of Ortiz et al. normalized to BP2000. The 
data from 14 MeV to 20 MeV are combined into a single bin. 
The horizontal solid line shows the measured total flux, while 
the dotted band around this line indicates the energy cor- 
related uncertainty. Error bars show statistical and energy- 
uncorrelated errors added in quadrature. 



in the day- night asymmetry measurement. The largest 
sources of error in the asymmetry are the energy scale and 
resolution (lo.oii) an< ^ the- non-flat background shape of 
the cos6* sun distribution (±0.004). 

Fig. shows the measured recoil electron energy spec- 
trum relative to the Ortiz et al. spectrum normalized to 
BP2000. (See Table |ffl] for numerical values.) A fit to an 
undistorted energy spectrum gives x 2 /d.o./. = 19.0/18. 
Energy-correlated systematic errors are considered in the 
definition of \ 2 ftoj- The energy-correlated systematic 



FIG. 4: Seasonal variation of the solar neutrino flux. The 
curve shows the expected seasonal variation of the flux in- 
troduced by the eccentricity of the Earth's orbit. Error bars 
show statistical errors only. 



error (shown in Fig. y as a band around the total flux) is 
due to uncertainties that could cause a systematic shift in 
the energy spectrum. The sources of this error are uncer- 
tainties in the energy scale, resolution, and the reference 
8 B spectrum against which the data are compared. 

The seasonal dependence of the solar neutrino flux 
is shown in Fig. [|. (See Table IV for numerical val- 
ues.) The points represent the measured flux, and the 
curve shows the expected variation due to the orbital 
eccentricity of the Earth (assuming no neutrino oscilla- 
tions, and normalized to the measured total flux). The 
data are consistent with the expected annual variation 
{x 2 /d.o.f. = 3.7/7). Systematic errors are included in 
the calculation of x 2 . The total systematic error on the 
relative flux values in each seasonal bin is ±1.3%, the 
largest sources coming from energy scale and resolution 
(l~i'i<y ) an d reduction cut efficiency (±0.6%). 

The hep neutrino flux given by BP2000 is 9.3 x 
10 3 cm~ 2 s _1 H, which is three orders of magnitude 



smaller than the 8 B neutrino flux. Although the theo- 
retically calculated hep flux is highly uncertain because 
of many delicate cancellations in calculating the astro- 
physical S factor, a recent calculation by Marcucci, et 
al. |2lj shows that the flux is unlikely to be drastically 
larger than the value given in BP2000. The effect of hep 
neutrinos on solar neutrino measurements at SK is ex- 
pected to be small. However, since the end-point of the 
hep neutrino spectrum is 18.77 MeV compared to about 
16 MeV for the 8 B spectrum, the high energy end of the 
-^recoil spectrum should be relatively enriched with hep 
neutrinos. An unexpectedly large hep flux may distort 
the -Erocoii spectrum. In our measurement of the hep flux, 
we extract the number of events in the window -Brecon = 
18-21 MeV from a cos sun distribution like the one shown 
in Fig. [l| This window was chosen because it optimizes 
the significance of the hep flux measurement in MC as- 
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FIG. 5: Energy spectrum of recoil electrons produced by B 
and hep neutrinos, in 1 MeV bins. The points show data with 
statistical error bars. The curves show expected spectra with 
various hep contributions to the best-fit 8 B spectrum. The 
solid, dotted, and dashed curves show the spectrum with 1, 
4.3, and times the BP2000 hep flux, respectively. 

suming BP2000 8 B and hep fluxes. We find 1.3 ± 2.0 
events in the chosen window. Assuming that all of these 
events are due to hep neutrinos, the 90% confidence level 
upper limit of the hep neutrino flux is 40 x 10 3 cm -2 s^ 1 
(4.3 times the BP2000 prediction). Fig. [H] shows the ex- 
pected energy spectra with various hep contributions. 
In summary, SK has lowered the analysis energy 



threshold to 5.0 MeV, collected more than twice the 
data previously reported, and reduced systematic er- 
rors through refinements in data analysis and exten- 
sive detector calibrations. With those improvements, 
and with the 18,464 observed solar neutrino events, SK 
provides very precise measurements of the recoil elec- 
tron energy spectrum, day-night flux asymmetry, and 
the absolute solar neutrino flux. The measured flux is 
45.1 ± 0.5 (stat.) tiA (sys.)% of the BP2000 prediction. 
We found no statistically significant energy spectrum dis- 
tortion (x 2 /d.o.f. = 19.0/18 relative to the predicted 8 B 
spectrum), and the day-night flux difference of 3.3% of 
the average flux is 1.3 er from zero. However, the precision 
of these measurements should provide strong and impor- 
tant constraints on the neutrino oscillation parameters. 
The seasonal dependence of the flux shows the expected 
7% annal variation due to the eccentricity of the Earth's 
orbit. This is the first neutrino-based observation of the 
Earth's orbital eccentricity A stringent limit on the hep 
neutrino flux ($hep < 40 x 10 3 cm~ 2 s _1 ) was obtained, 
which corresponds to 4.3 times the predicted value from 
BP2000. 
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Day-Night 


Zenith angle 


Data/SSM 


Si 


DAY 


-1 < cos0 z < 


u -^°-0.007 






MANl 
MAN2 
MAN3 
MAN4 
MAN5 
CORE 


< cos 0* < 0.16 
0.16 < cos0 z < 0.33 
0.33 < cos0 z < 0.50 
0.50 < cos0 z < 0.67 
0.67 < cos0 z < 0.84 

0.84 < cos0 z < 1 


n 47C+0.020 

U-*'O_0 .020 

47fi+ - 018 

476+ ' 016 
0.447lg : gl? 
0.438lHi 9 9 


+1.3 % 
-1.2 /0 
+ 1.3 % 
-1.2 /0 
+ 1-30/ 
-1.2 /o 
+ 1.3q/ 
-1.2 /o 
+ 1.30/ 

+ 1.307 

-1.2/0 



TABLE II: Numerical results of the day/night analysis. The zenith angle region (2nd column), the ratio of observed and 
expected number of events (3rd column), and la error of systematic error (4th column). Systematic errors are relative to 
Day-flux. Energy range is 5.0-20.0MeV. 



Energy 


Data/SSM 


Si, cor. 


$i,uncor. 


(MeV) 


ALL 


DAY 


NIGHT 






5.0-5.5 





* oo -0.046 


0.438t°;°» 





404+0. 063 
^°^-0.063 


+ O.207 
-0.2/0 


+3 
-3 


\% 


5.5-6.0 





,,0.0+0.024 
4 ^ S -0.024 


0.4281°;°* 





446 +0.034 
^"-0.034 


+ O.207 
-0.2/0 


+1 
-1 


1% 


6.0-6.5 





,or+0.019 
^°°-0.019 


426+ - 027 





444+0.027 

*^*-0.027 


+ 0.3O/ 
-O.3/0 


+1 
-1 


1% 


6.5-7.0 





400+0.015 
4^8_o.oi5 


431+0 021 

O.^01_o.021 





444+0.021 

^* -0.021 


+0.5O/ 
-0.6 /o 


+1 
-1 


t% 


7.0-7.5 





4 63I .oi5 


462+ 022 





464 +0.022 
* oy: -0.022 


+O.807 
-0.8 /o 


+1 
-1 


1% 


7.5-8.0 





,,0 0+0.016 
4S 3_ .016 


4 94+ - 023 

U.1»'±_ 023 





4^9+0.022 
^' z -0.022 


+ 1.0 % 
-1.1 /0 


+1 
-1 


1% 


8.0-8.5 





4fit -+0.017 


O.4:0^_o.o23 





477+0.023 
' -0.023 


+ 1.40/ 
-I.3/0 


+1 
-1 


1% 


8.5-9.0 





^00+0.017 
4^»_o.oi7 


402+ - 024 

U. i iUZ_ 024 





470+O.O24 
4: '°-0.024 


+ 1.707 
-I.7/0 


+1 
-1 


1% 


9.0-9.5 





4C.A+0.018 
4dU -0.018 


u.tot_ 026 





446 +0.025 
^"-0.025 


+2.I07 
-2.0 /o 


+1 
-1 


1% 


9.5-10.0 





455+0.019 
^ dd -0.019 


44Q+ - 027 
u.fty_ 027 





4fif) +0.027 
*«O_ 027 


+2.50/ 
-2.3/0 


+1 
-1 


1% 


10.0-10.5 





442 +o.02i 

44Z -0.021 


0.430tg : °» 





454+0.029 

^ J *-0.029 


+3.00/ 
-2.7 /0 


+1 
-1 


\% 


10.5-11.0 





4H7+0.022 
+ u ' -0.022 


0.386iro 3 3°o 





426 +0.032 
4:ZU -0.032 


+3.4o/ 
-3.2/0 


+1 
-1 


1% 


11.0-11.5 





4rr+0.026 
^ dd -0.026 


439+ 036 





470 +0.037 
^' u -0.037 


+3.9 % 
-3.6 /0 


+1 
-1 


\% 


11.5-12.0 





490+0. 028 
^ ZO -0.028 


O.^OO_ 042 





QQ4+0.038 


+4.50/ 


+1 
-1 


1% 


12.0-12.5 





499+0.033 
^ zz -0.033 


0.3891°;°* 





4CC+0.047 
4:dd -0.047 


+5.10/ 
-4.8 /0 


+1 
-1 


\% 


12.5-13.0 





401 +0.041 
481_ 041 


0.514tro 6 6 J 





451+0.055 

4di -0.055 


+5.807 
-5.4/0 


+1 
-1 


1% 


13.0-13.5 





401+0.047 

^ O1 -0.047 


0.468_ 070 





OQ7+0.063 
aS ' -0.063 


+6.50/ 
-6.2/0 


+1 
-1 


1% 


13.5-14.0 





Rn o+0.065 
«OO_ 065 


n m +0.092 
U.OOl_ 092 





ccq+0.094 
uda -0.094 


+7.4o/ 
-7.0 /o 


+1 
-1 


1% 


14.0-20.0 





4QO+0.049 
^ yo -0.049 


4 qo+ 067 





550+0.071 

ddy -0.071 


+ 10.7(V 
-9.5 /o 


+1 
-1 


1% 



TABLE III: Numerical results of the energy spectrum analysis. The ratio of observed and expected number of events in all-time 
(2nd column), in day-time (3rd column), in night-time (4th column), la error of correlated systematic error (5th column), and 
la error of uncorrelated systematic error (6th column). Systematic errors are relative. 
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Time period 


Data/SSM 


Data/SSM at 1 AU 




Jan. 1 - Feb. 16 
Feb. 17 - Apr. 2 
Apr. 3 - May 19 
May 20 - Jul. 4 
Jul. 5 - Aug. 19 
Aug. 20 - Oct. 4 
Oct. 5 - Nov. 18 
Nov. 19 - Dec. 31 


0.475t°:°S 
0.456lHi 4 4 

o.44it° : °i° 

0.438t° : °lS 
426+ 013 
462+ 013 
0.439+g-gii 
0.458lg : gi« 


o.46ii :° s 
o.45oi2:2i4 

446+ ' 015 
n 4ki +0-014 

U.4Ol_ 014 
4QQ+0 013 

46Q+ 014 

<J-^O»_0 .014 


+ 1.3o/ 
-1.3 /o 
+ 1.3o/ 
-1.3 /o 

-1.3 /o 
+ 1.3q/ 
-1.3 /o 
+ 1.30/ 
-1.3 /o 
+ 1.3 % 
-1.3 /0 
+ 1.3 % 
-1.3 /0 
+ 1.3 % 
-1.3 /0 



TABLE IV: Numerical results of the seasonal analysis. The time period (1st column), the ratio of observed and expected 
number of events without eccentricity correction (2nd column), the ratio of observed and expected number of events with 
eccentricity correction (3rd column), and la error of systematic error (4th column). Systematic errors are relative. Energy 
range is 5.0-20.0MeV. 



